Glissile dislocations with transient cores in silicon 
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We report an unexpected characteristic of dislocation cores in silicon. Using first-principles cal- 
culations, we show that all the stable core configurations for a non-dissociated 60° dislocation are 
sessile. The only glissile configuration, previously obtained by nucleation from surfaces, surprins- 
ingly corresponds to an unstable core. As a result, the 60° dislocation motion is solely driven by 
stress, with no thermal activation. We predict that this original feature could be relevant in situ- 
ations for which large stresses occur, such as mechanical deformation at room temperature. Our 
work also suggests that post-mortem observations of stable dislocations could be misleading, and 
that mobile unstable dislocation cores should be taken into account in theoretical investigations. 
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Dislocations are linear defects present in most of mate- 
rials, and have been largely studied since they are known 
to strongly influence many properties, primarily mechan- 
ical and electrical [H, 0] ■ Dislocations often play an im- 
portant role in the plastic deformation of materials, or in 
epitaxial growth of thin films, because they allow to effi- 
ciently relax stress when they move. Mobility is therefore 
a key property of dislocations, in connection with macro- 
scopic characteristics such as brittleness or ductility. In 
a lattice, at rest, dislocations lie in Peierls valleys sepa- 
rated by energy barriers, whose heights define the lattice 
resistance to dislocation motion 3]. In materials with 
deep Peierls valleys, i.e. with moderate to large lattice 
resistance, such as bcc metals and semiconductors, the 
properties of dislocations and in particular their mobil- 
ity typically depend on their core structure, which can 
be complex and reconstructed In order to move, a 
dislocation must overcome the energy barriers, what can 
be achieved thanks to the combined action of stress and 
temperature. 

Considering first the role of temperature, the motion 
of dislocations is occurring by thermally activated pro- 
cesses such as creation and migration of kinks or jogs. 
With temperature alone, all possible directions for dis- 
placement are equivalent, and a moving dislocation would 
behave like a random walker. Conversely, in presence of 
an applied stress, dislocation motion along a specific di- 
rection will be favored. In addition, stress will reduce 
the energy barriers, thus making easier the dislocation 
displacement. An important quantity related to dislo- 
cation mobility is the Peierls stress, which is the mini- 
mum stress required to move a straight dislocation along 
a given direction without any thermal activation. Dis- 
locations moving in the usual range of temperature and 
applied stress are called glissile, whereas non-moving dis- 
locations are called sessile. 

Theoretical characterizations of dislocation mobility 
are usually done according to the following procedure. 
First, low energy configurations for a dislocation core are 



determined. Then, for the most stable ones, the Peierls 
stress and the role of thermal activation are investigated. 
Such an approach is based on the hypothesis that the core 
configuration for a dislocation at rest is the same than 
for a moving dislocation. This is a reasonable and gen- 
eral assumption, also employed for point defects diffusion 
for instance, that has always been verified as far as we 
know. However, in this letter, we reveal a contradictory 
situation, where all stable configurations for a dislocation 
core are sessile, whereas the only glissile configuration is 
unstable, i.e. not lying in a Peierls valley. Such an un- 
expected result is fully contrasting with all previous the- 
oretical investigations of dislocation cores, and suggests 
that it could be advisory to determine mobile dislocation 
core rather than stable ones. This result also hints that 
the transient character of some dislocation cores would 
prevent a direct observation with usual techniques such 
as post-mortem transmission electron microscopy. 

In this work, we have used silicon as a model of ma- 
terial with high lattice friction. In its cubic diamond 
structure, usual dislocations have Burgers vector b equal 
to ao/2(110), ao beingthe lattice parameter, and orien- 
tations screw and 60° @]. At high temperature (> 700K) 
and low stress, both are dissociated into 30° or 90° par- 
tial dislocations, whereas at low temperature (< 700K) 
and high stress (~ 1.5 GPa) they do not dissociate 0,0). 
In the latter regime, compared to the screw, the 60° dis- 
location is commonly assumed to be more mobile. As 
such, it is expected to play an important role for relax- 
ing epitaxial stresses in thin films 0- Also, it has been 
shown that the 60° dislocation can form from surface de- 
fects [H H, [Toll . a sharp corner [ll[ , or in the vicinity of 
crack fronts [12| , highlighting the fundamental role of this 
dislocation in the plastic deformation of silicon-based ma- 
terials. The actual knowledge of the 60° dislocation core 
geometry is mainly based on an early structure analy- 
sis, suggesting two possible core configurations. The first 
one, proposed by Hornstra [lH and called SI in this work 
(Fig. [T]), is located in shuffle {111} planes. It is charac- 
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TABLE I: Energy differences, for possible core configurations 
of a 60° dislocation. The most stable reconstructed core con- 
figuration G is taken as reference. 
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FIG. 1: Possible atomistic configurations for a 60° dislocation 
core in silicon, following the notation described in the text. 
Bonds are drawn according to a criterion of distance between 
atoms. In the case of the S2 configuration, the dashed line 
indicates a distance of 2.82 A between two atoms in the core. 
An ellipse allows to mark out the two atoms with maximum 
displacements during the SI— >S3 transformation. 



terized by a 8-atoms ring, containing an atom carrying 
a dangling bond. Classical potential calculations showed 
that this dislocation core is stable and that the Peierls 
stress is about 1.2 GPa 14]. This SI configuration has 



been identified in previous investigations of the nucle- 
ation of dislocations in silicon [8l 



lfj, [ll| • It has also 



been selected for studying the interaction between dislo- 
cation and intrinsic point defects [IH . A second possible 
configuration (called G in Fig.[T|) is located in glide {111} 
planes 0]. In that case, the core structure is character- 
ized by two rings of 5 and 7 atoms, and the absence of 
dangling bonds. Although more stable than the SI core, 
this configuration has been shown to be sessile when con- 
sidered for relaxing epitaxial stresses in thin films [7j. 

Using state of the art electronic structure calculations, 
we have investigated the properties of 60° dislocation, in 
order to ensure earlier results on the stability and mo- 
bility of this dislocation. In a first step, we generated 
undissociated 60° dislocation core configurations in a sil- 
icon cubic diamond crystal using displacements according 
to anisotropic elasticity theory [16[ . The dislocation cen- 
ter was set to various positions relatively to the lattice, in 
order to produce several different dislocation cores. All 
generated configurations were then relaxed using a clas- 
sical interatomic potential Only two stable disloca- 
tion core structures were finally obtained corresponding 





SI 


S2 


S3 


AE (eV/6) 


~2.4 


1.4 


0.6 



to the configurations G and SI, shown in the fig.[T] Both 
structures are in agreement with previous works, showing 
similar topological features as previously depicted. 

Next, we investigated the relative stability of both con- 
figurations by performing total energy calculations in the 
framework of density functional theory. Initial compu- 
tations were made using the DFTB+ code, relying on 
the tight-bin ding approximation and an appropri- 
ate basis set 1191 . We then employed the first principles 
SIESTA code [20( for obtaining highly accurate results. 
In the latter case, the local density approximation, norm- 
conserving pseudoptentials [21|, and double-^ polarized 
localized orbitals basis set with a cutoff of 6 bohr were 
used. Optimized lattice constants of 5.46 A for DFTB+ 
and 5.404 A for SIESTA were obtained, in good agree- 
ment with the experimental value [22I ]. Dislocations in 
bulk silicon were mainly modelled using a quadrupolar 
arrangement together with periodic boundary conditions, 
the specific geom etry of the cells including only two dislo- 
cation cores [23[. Cell dimensions were 12 x 12 along [121] 
and [111] axis, ensuring that the two dislocation cores 
were separated by 6 hexagons along [121], and 1 or 2 pe- 
riods along the dislocation line orientation [101], leading 
to calculations with either 144 or 288 atoms. Depending 
on the system size, 1 or 2 k-points along the dislocation 
line were used. In order to check the possible effect of 
boundary conditions, test calculations with DFTB+ and 
clusters including a single dislocation core and at most 
232 Si atoms were also performed, surface atoms being 
fixed and saturated with hydrogen. Both DFTB+ and 
SIESTA results were fully in agreement, that is why only 
SIESTA results are described below. 

Starting from an initial configuration G, wc found very 
little structural modifications after forces relaxation, per- 
formed until all forces were below 10 -2 eV/A. However, 
we found that the energy of this specific configuration 
could be further lowered by 0.7 eV/6 if the core is recon- 
structed with a double period along the dislocation line. 
This reconstructed core has been previously proposed in 
the case of 60° dislocation in diamond [24]. In the case of 
an initial shuffle SI core, we found the surprising result 
that this geometry is not stable and transforms to an- 
other configuration, approximately 1 eV/b lower, which 
we call here S2 (fig. Q}. This unstable behavior, non oc- 
curing with classical potentials [3], was obtained in all 
electronic structure calculations, performed either with 
DFTB+ and SIESTA. It also did not depend on bound- 
ary conditions, since simulations done with cluster sys- 
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terns lead to a similar outcome. We found that it is pos- 
sible to further decrease the dislocation core energy by 
slight atomic rearrangements. The new core structure, 
called S3, is shown in the fig.[TJ and lowers the energy by 

0. 8 eV/b compared to the S2 configuration. Compared 
to the glide core G with the double period reconstruction 
along the dislocation line, the S3 configuration is 0.6 eV/6 
higher in energy. The description of the geometry of new 
configurations S2 and S3, as well as of the transformation 
mechanisms, is not the focus of this letter, and will be 
reported elsewhere. Summarizing our results in Table U 
we found two original stable configurations S2 and S3 
for shuffle cores, whereas the usually considered geome- 
try SI is shown to be unstable. The problematics of 60° 
dislocation in silicon is then more complicated than first 
thought. 

As stated previously, an important property of a dislo- 
cation is its mobility, which can be characterized by the 
Peierls stress. We performed calculations of the Peierls 
stress considering all possible stable 60° dislocation con- 
figurations. An increasing shear stress, either parallel 
(defined as positive) or anti-parallel (negative) to the 
Burgers vector, is applied by straining the computational 
cell, keeping a constant volume cell, until forces relax- 
ation leads to a displacement of the dislocation. Disloca- 
tion cores not displaced for shear strains larger than 15%, 

1. e. close to the ideal shear strain 25, 2(1, are considered 
sessile. Note that parallel and anti-parallel shear stress 
orientations are not equivalent for a 60° dislocation in 
the cubic diamond structure, yielding two values for the 
Peierls stress. First, for shear strains as large as 15% and 
for both orientations, we found that there is no displace- 
ment of the glide core G. This configuration is therefore 
sessile, in agreement with previous investigations 0. No 
attempts were made for determining the shear stress of 
the double period glide core, but it is very likely that this 
highly stable and reconstructed core is even more diffi- 
cult to move. Then, we considered both S2 and S3 shuffle 
cores. The former is transformed to a S3 configuration 
for an applied shear strain of -2% (1.1 GPa), whereas it 
remains still for all positive strains. The S3 configuration 
is also not displaced over all the strain range, which is 
not surprising since it corresponds to a very stable and 
highly reconstructed core. Additional finite temperature 
calculations, to be reported elsewhere, indicated that the 
thermally activated motion of both G and S3 dislocation 
cores was not favored. 

In previous works, it was commonly assumed that 60° 
dislocations could exist in two stable configurations, G 
and SI, with only the latter being glissile. Our sim- 
ulations rather indicate that the shuffle core SI is un- 
stable, and that it can be transformed in two different 
structures, S2 and S3, both being sessile. As a matter 
of fact, we found no stable and glissile 60° dislocation 
cores. It is a surprising result, because it is difficult to 
imagine dislocation loops in silicon with no mobile 60° 



segments. Furthermore, previous investigations of dislo- 
cation nucleation from silicon surfaces, made with first 
principles calculations, point at formation and propaga- 
tion of 60° dislocations, with core structures seemingly 
close to a SI configuration Q. This last result raises 
the question of the existence of dislocation cores which 
could be mobile albeit unstable. To tackle this issue, 
we tried to determine the stress required for displacing 
the unstable SI core in the silicon lattice. We used an 
initial SI core obtained as a stable configuration in clas- 
sical potential calculations, which was subsequently first 
strained, then relaxed with first principles simulations. 
For applied strains lower than 2%, we found that the SI 
core evolves to the S2 or S3 configurations. However, for 
larger strains, typically between 2% and 5% [53] (1.1 to 
2.8 GPa), the SI core is displaced until it encounters the 
other dislocation in the cell. During its displacement, the 
geometry of the SI core is approximately preserved. Ad- 
ditional calculations have been performed using cluster 
cells and DFTB+, leading to similar results. This re- 
sult is in agreement with previous investigations . Our 
conclusion is that only an unstable 60° dislocation core 
is mobile. This glissile core could therefore be described 
as transient, since it exists only in motion. 

From our calculations, we propose that 60° disloca- 
tions that are assumed to participate to the deformation 
of bulk semiconductors in the high stress/low tempera- 
ture regime, or to the relaxation of large stresses in epi- 
taxial films, or to the nucleation of dislocations from sur- 
faces, are characterized by this transient core. Typically, 
there is first nucleation of a 60° dislocation in the unsta- 
ble mobile configuration followed by the immediate glide 
of the dislocation, if the applied stress is larger than the 
Peierls stress. If the stress decreases below this thresh- 
old, the dislocation will stop moving and will relax to the 
sessile core configuration. Such a scenario would occur 
for any reasons that put the dislocation motionless, even 
during a very short time. This result, in agreement with 
previous works showing the nucleation and propagation 
of the 60° dislocation in silicon Q , has several important 
and unusual implications. First, the stress required to 
move a 60° dislocation should always be larger than the 
Peierls stress. As a consequence, there is no thermally 
activated motion of the 60° dislocation, even at finite 
temperature. Second, the role of dislocation nucleation 
is essential, since it controls the whole process. Finally, 
since 60° dislocations coming to a halt become sessile, 
new dislocations have to be formed from other sources 
for relaxing excess or additional stresses. 

In a more general perspective, here we show that 
mobile dislocation cores, allowing to relax mechanical 
stresses, can be intrinsically unstable. Consequently, 
they can hardly be observed experimentally, and espe- 
cially not in ex-situ experiments. It also raises issues 
regarding the theoretical investigation of dislocation sta- 
bility and mobility. In fact, usually one first determines 
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the stable configuration of the dislocation in the Peierls 
valleys, then after how this configuration would move as 
a function of stress and temperature. Instead, in our 
case, dislocation core mobility, and not stability, should 
be investigated first. Note that it could be a difficult or 
even unfeasible task. Here, it was fortunate that classi- 
cal potential calculations spuriously allow to stabilize the 
mobile core. 



Finally, it is interesting to discuss if this result is a 
rarity, specific to silicon, or whether it could occur in 
other systems. Certainly, an important condition is that 
several core configurations should exist for a given dis- 
location. This is certainly true in materials where re- 
construction of the dislocation core could occur, leading 
to several more or less complex and stable configurations. 
All materials showing a covalent bonding character could 
therefore be concerned: usual semiconductors (Group IV, 
III-V, or II- VI) or more complex systems such as min- 
erals present in the Earth mantle, or ceramics. Also, 
stress should be the main cause for dislocation motion 
compared to thermal activation, suggesting that tran- 
sient glissile cores could occur in the high stress / low 
temperature regime. Plastic deformation of several com- 
pound semiconductors in this regime has been shown to 
be similar to what is known for silicon 28, l29|. Tran- 



sient glissile dislocations could also be of prime impor- 
tance for two important domains of material sciences. 
The first concerns the onset of plasticity in materials with 
one or more nanometric dimensions. In such systems, in 
which there are often no pre-existing defects, and disloca- 
tion multiplication is hindered, a large amount of stress 
can be stored in thin films, nanopillars or nanowires dur- 
ing mechanical testing at room temperature. After nu- 
cleation of a primary dislocation from surfaces, disloca- 
tion avalanches lead to sudden and large relaxation of 
stress. It is interesting to note that these stresses are 
typically a significant fraction of the ideal strength of 
the bulk material, thus large enough for transient glissile 
core to play a role. As reported previously, such an oc- 
curence has been shown for a silicon thin film under large 
stress Q. We believe that this scenario could also occur 
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30l 31, 32, 33j. The second situation con- 
cerns brittle-to-ductile transition, for which dislocation 
nucleation in the vicinity of the crack front is a critical 
factor. Here again, large stresses are present, and the 
occurence of transient glissile dislocation cores should be 
considered. 
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